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Abstract 

In this paper, we investigate the vahdity of the generalized second law (GSL) of thermodynamics 
in flat FRW chameleon cosmology where the boundary of the universe is assumed to be enclosed by 
the dynamical apparent horizon. It has been shown that, in a bouncing scenario for the universe 
with phantom crossing, the total entropy decreases with time in the contracting epoch, whereas, 
the dynamics of the internal and horizon entropies depends on the behavior of both equation of 
state and bubble parameters. 
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1. INTRODUCTION 



Motivated by the black hole physics, it was realized that there is a profound connection 
between gravity and thermodynamics. In Einstein gravity, the evidence of this connec- 
tion was first discovered in [l| by deriving the Einstein equation from the proportionality 
of entropy and horizon area together with the first law of thermodynamics in the Rindler 
spacetime. For a general static spherically symmetric spacetime, Padmanabhan pointed out 
that Einstein equations at the horizon give rise to the first law of thermodynamics j^. Re- 
cently, the study on the connection between gravity and thermodynamics has been extended 
to cosmological context. Frolov and Kofman jsl employed the approach proposed by Jacob- 
son jl| to a quasi-de Sitter geometry of inflationary universe, and calculated the energy flux 
of a background slow-roll scalar through the quasi-de Sitter apparent horizon. By applying 
the first law of thermodynamics to a cosmological horizon, Danielsson obtained Friedmann 
equation in the expanding universe Q. In the quintessence dominated aecelerating universe, 
Bousso [5| showed that the first law of thermodynamics holds at the apparent horizon. The 



relation between gravity and thermodynamics 



itationa 
gravity 



. theories, including Lovelock gravity 



las been further disclosed in extended grav- 
braneworld gravity 8|-j9|, nonlinear 



etc. In scalar-tensor gravity, it was ar- 



12| and investigating 



13| have also been the 



7|, [lO|-[ll|, scalar-tensor gravity [7|- 
gued that the non-equilibrium thermodynamics instead of the equilibrium thermodynamics 
should be taken into account to build the relation to gravity 10|-ll|. Determining ther- 
modynamic parameters of an (accelerated) expanding universe and verifying the first and 
the second thermodynamics law for different cosmological horizons 
the relation between dynamics and thermodynamics of the universe 
subjects of many researches in recent years. In particular, the validity of the generalized 
second law (GSL) [l^ which state that entropy of the fluid inside the horizon plus the en- 
tropy associated with the apparent horizon do not decrease with time, has been the subject 
of many studies. In this paper, we study the thermodynamics properties of the chameleon 
cosmology in which the universe undergoes a kind of bouncing phase transition and phantom 
crossing. In section two, we drive the field equations for the model. The thermodynamics 
properties of the model is also studied in section three. 
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2. THE MODEL 



In this section, we consider the chameleon gravity in the presence of matter with the 



action given by 



15| 



S = I [-^ - l^^" + V{<jy) + /(0)L„]v^dx^ (1) 



where R is Ricci scalar, G is the newtonian constant gravity and is the chameleon scalar 
field with the potential V{(f)). Unlike the usual Einstein-Hilbert action, the matter La- 
grangian Lm is modified as f{(f))Lm, where /(</>) is an analytic function of (p. This last term 
in Lagrangian brings about the nonminimal interaction between matter and chameleon field. 
The variation of the action ([T]) with respect to the metric tensor components in a spatially 
fiat FRW cosmology yields the field equations, 

= p„/ + ^02 + ^(0), (2) 

2H + 3H^ = -jp^f - + Vi<P), (3) 

where we put SnG = c = h = 1 and assume a perfect fiuid with = ■ypm- The dots means 
derivatives with respect to the cosmic time t. The energy density pm is the matter energy 
density in the universe. Also variation of the action ([T]) with respect to the scalar field 
provides the wave equation for the chameleon field as, 

0(0 + 3/f0) = - 1(1 - 3^)pJ. (4) 

From equations ([2]), ([3]) and (jl]), one can easily arrive at the relation, 

ipj) + 3H{1 + 7)p„/ = lp„(l - 37)0/, (5) 

which readily integrates to yield 

- ^ 

with M as a constant of integration. From equations (|2]) and ([3]) and in comparison with 
the standard friedmann equations we identify pe// and Pefj as, 

Peff^Pm.f+l<P' + Vi<l>). (7) 
Peff^lPmf+l^^-V{(P), (8) 



with the equation of state, peff = UejfPeff- From equations ([7]) and ([8]), one leads to, 
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-Y[P™/(7 - ^eff) - V{1 + Ueff)]. 



(9) 



In the next section we investigate the thermodynamics and GSL of the modeL 



GSL 



In the following we make two assumptions: 1) In addition to the entropy of the universe 
inside the horizon, an entropy is associated to the apparent/event horizon. 2) with the local 
equilibrium hypothesis, the energy would not spontaneously flow between the horizon and 
the fluid inside the horizon, the latter would be at variance with the FRW geometry. 

According to the GSL in an expanding universe, entropy of the viscous DE, DM and 
radiation inside the horizon plus the one associated with the apparent horizon do not decrease 
with time. However, in the cosmological models where the universe has a bouncing behavior, 
the rate of change of the entropy can be negative as will be discussed shortly. 

In general, there are two approaches to validate the GSL on apparent /event horizons: i) 



by using first law of thermodynamics and find entropy relation on the horizons i.e. 



ThdSh = -dEh = ^TiRlHT^.K^'f^^dt = A7iRlH{p,ff+p,ff)dt, 



(10) 



where 



Hr, 0, 0) are the (approximate) Killing vector ( the generators of the hori- 



zon), or the future directed ingoing null vector field 



and "h" stands for the horizon, ii) 



in the field equations, by employing the horizon entropy and temperature formula on the 



horizon 
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Sh 



ttR 



(11) 
(12) 



Note that only on the apparent horizon the two approaches are equivalent 19|. Furthermore, 
the recent observational data from type la Supernovae suggests that in an accelerating 
universe the enveloping surface should be the apparent horizon rather than the event one 



20|. Therefore, from now on, we assume that the universe is enclosed by the dynamical 



apparent horizon with the radius given by Rh 



2l| 
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By the horizon entropy and temperature given in equations (fTT!) and (fT2|) . the dynamics 
of the entropy on the apparent horizon is ij] , 

Sh = 27iRhRh. (13) 



Also, from the Gibbs equation, the entropy of the universe inside t 



le horizon can be related 



to its effective energy density and pressure in the horizon with 22| 



TdS,n=PeffdV + d{Ei^), (14) 

where Sin is the internal entropy within the apparent horizon and Pe// is the effective pressure 
in the model. If there is no energy exchange between outside and inside of the apparent 
horizon, thermal equilibrium realizes that T = T^. Here, the expression for internal energy 
can be written as Ein = peffV, with V = |vri?|. Therefore, from equation f|T^ . by using 
Friedmann equation and after doing some calculation we find that the rate of change of the 
internal entropy, horizon entropy and total entropy are respectively. 

Sir, = 127rX^(l + ^e//)(l + 3a;e//), (15) 

Sh = 247T^RlH{l+Ueff), (16) 

St = 367r^RlH{l + uj,ffy. (17) 

As already noted the GSL states that the total entropy is not a decreasing function of time, 
or St > 0. However, if we assume a bouncing universe as in our chameleonic model, then we 
expect a decreasing total entropy when the universe is in contracting epoch. 

For the rate of change of the internal entropy, equation ( fTSl) . we find that when H < 
and Ueff < —1 or Uefj > —1/3, then Sin < 0. Also, when H > and —1 < Uefj < —1/3, 
again Sin < 0. On the other hand, when H < and — 1 < We// < —1/3 or H > and 
Ueff < — 1 or Ueff > —1/3, we obtain that Sin > 0. we conclude that the behavior of the 
internal entropy depends both on the bubble parameter and effective EoS parameter of the 
model. From equation f|T6l) . it can be seen that when H and (1 + ooeff) have the same sign, 
then Sh > 0. Otherwise, Sh < 0. Finally, in equation ( |T71) . the sign of the total rate of 
change of the entropy, St, is independent of the We//, while depending only on the sign of 
the bubble parameter. For an expanding/contracting universe, we have St ^ 0. 

In the following we will numerically investigate the validity of GSL on Chameleon model. 
In order to close the system of equations we make the following ansatz: We take both 



potential and / behave exponentially as V{(f>) = Voe^^'^ and /(0) = foe^^"^ where Si,S2,Vo 
and /o are arbitrary constants. There are no priori physical motivation for these choices, so 
it is only purely phenomenological which leads to the desired behavior of phantom crossing. 
The parameters 6i and 62 are dimensionless constants characterizing the slope of potential 
V{(f)) and Here, we assume that the universe is filled with cold dark matter, i.e. 7 = 

The graphs in Fig.l provide a dynamical universe with contraction for t < 0, bouncing 
at t = and then expansion for t > 0. From equation (|9]), one can see that for 7 = the 
crossing occurs only for negative /. Also we see that at the bouncing point where the scale 
factor a{t) is not zero we avoid singularity faced in the standard cosmology 15|. In Fig. la), 
we observe that the effective EoS parameter, Ueff crosses the phantom divide line twice at 
t = 0.75 and t = 1.25. 




.75 -0.5 -0.25 0.25 0.5 0.75 I 1.25 1.5 1.75 2 





Fig. 1: The behavior of the effective EoS parameter, We//? hubble parameter, H{t), 
and the scale factor a(t) with respect to the cosmological time t. The potential V and / are 
V{(f)) = Vocxp and /(0) = /oexp (^20) where /o = —10, Vq = 15 and 6i = —1, 62 = —I. 

I.e. 00 = 1, 0(0) = -0.2. 

In Fig. 2) the rate of change of the internal entropy, horizon entropy and total entropy 
with respect to the effective EoS parameter in our model are shown. We also show the 
dynamics of the H{t) for a comparison with the analytical findings. For example, for the 
internal entropy, from the graph, and for H < 0, one observes that for Ueff > —1/3, then 
Sin < 0; for —0.6 < u^ff < —0.3, then Sm > 0; and for —1 < u^ff < —0.6 when H > 
again Sin < 0. Similar argument can be made for horizon and total entropies and as a result 
it shows that the numerical calculation is in complete agreement with the analytical findings. 
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Fig. 2: The dynamics of the Sm, Sh, and St with respect to the effective EoS parameter, We//- 

In Fig. 3, for Wg// < —1, we see that Sh < 0, whereas, both Sin and St are positive. The 
negative Sh is due to the phantom behavior of the dark energy. At Wg// = —1, as can be seen 
all three S^, Sh, and St vanish, whereas, at Wg// = —1/3, only S^ vanishes, as expected. 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



Fig. 3: The dynamics of Sin, Sh and St with respect to the cosmological time 
and in comparison to the dynamics of the effective EoS parameter, cug//, 

Fig. 4 shows that for 7 = the crossing occurs only when 0^ and — p„/ are equal and at 
t = tcross = 0.75, 1.25, St = 0. Also when H = a.t t = 0, again St vanishes. 
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4. CONCLUSION 



In this paper, we study the thermodynamics of the chameleon cosmological model in 
which a light scalar field (chameleon field) nonminimally coupled to the matter Lagrangian. 
We find that the evolution of the scale factor of the universe is non-singular in a bouncing 
cosmology, with an initial contracting phase which lasts until to a non-vanishing minimal 
radius is reached and then smoothly transits into an expanding phase which provides a pos- 
sible solution to the singularity problem of standard Big Bang cosmology. We also show 
that the the total entropy of the universe in the model with bouncing behavior and dy- 
namical effective EoS parameter increases/decreases with time in the expansion/contraction 
period. However, for the internal and horizon entropies, the sign of Sin and Sh depend on 
the dynamics of both hubble and effective EoS parameters in the model. Further, even in 
an expanding universe, when —1 < Wg// < —1/3, the internal entropy decreases with time 
while the apparent horizon and total entropies do not. Another result we obtained is that 
in an expanding universe in phantom era where We// < —1, the horizon entropy decreases 
with time which is due to the phantom behavior of the dark energy. Finally, we would like 
to mention that the thermodynamics of the bouncing universe with nonminimal coupling is 
an interesting topic. From a thermodynamical point of view, due to the modification of the 
theory of gravity in the early universe and its effects to the universe evolution, the subject 
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deserves more attention in both theoretical and observational aspects. 
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